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We study the effects of dynamical imperfections in quantum computers. By considering an explicit 
example, we identify different regimes ranging from the low-frequency case, where the imperfections 
can be considered as static but with renormalized parameters, to the high frequency ffuctuations, 
where the effects of imperfections are completely wiped out. We generalize our results by proving a 
theorem on the dynamical evolution of a system in the presence of dynamical perturbations. 

PACS numbers: 03.67.Lx, 05.45.Mt, 24.10.Cn, 03.67.Mn, 03.67. 



In any experimental implementation of a quantum in- 
formation protocol one has to face the presence of 
errors. The coupling of the quantum computer to the sur- 
rounding environment is responsible for decoherence 
which ultimately degrades the performances of quantum 
computation. The presence of static imperfections, al- 
though not leading to any decoherence, may be also detri- 
mental for quantum computers. For instance, a small 
inaccuracy in the coupling constants, inducing as a con- 
sequence to errors in quantum gates, can be tolerated 
only up to a certain threshold [2j. Moreover, the role 
of static imperfections depends on the regime, chaotic or 
not, of the system under consideration [2j. The stability 
of a quantum computation in the presence of static im- 
perfections has been already analyzed both in terms of 
fidelity H, ^ |^ and entanglement . 

A strict separation in "static" imperfections and "dy- 
namical" noise may not be always satisfactory. Dynam- 
ical noise may be considered at the same level as static 
imperfections, if its evolution occurs on a scale much 
larger than the computational time. In Ref. Q it was 
suggested that the effects of static imperfections can be 
more disruptive than noise for quantum computation. In 
this Letter, we intend to explore this problem in more 
details. The model we consider, in spite of its simplicity, 
enables one to grasp the interplay between the differ- 
ent time scales that appear in the problem. We consider 
each qubit coupled to a stochastic variable which changes 
in time with a fixed frequency. Below a given thresh- 
old (frequency), the errors can be considered as static, 
and thus can be corrected by using any of the known 
methods. The difference between the chaotic and the 
other dynamical regimes, found for static imperfections, 
holds also in the quasi-static case. We then generalize 
our results, by proving a theorem that states that, under 
general assumptions, in a perturbed system, unitary dy- 
namical errors are averaged to zero in probability. Our 
results might be relevant in the context of the strate- 
gies that have been proposed during the last few years in 
order to suppress decoherence . 

Model - Following 0, 0, we model a quantum com- 



puter as a lattice of interacting spins (qubits). Due to 
the unavoidable presence of imperfections, the spacing 
between the up and down states (external field) and the 
couplings between the qubits (exchange interactions) are 
both random and fluctuate in time. We consider n qubits 
on a two-dimensional lattice, described by the Hamilto- 
nian 

n 

H(t) = 5][Ao + <5,(i)]ap^ + J^,{i)a^ol^\ (1) 

where the CTq^'s (a — x,y,z) are the Pauli matrices for 
qubit i and the second sum runs over nearest-neighbor 
pairs. The energy spacing between the up and down 
states of a qubit is Aq -I- 5i (t) , where the 6i {t) 's are uni- 
formly distributed in the interval [—6/2,6/2] and the 
Jij{t)'s in the interval [—J, J] (zero means and vari- 
ances 6^a^ and 4J^ct^, respectively, with cr^ = 1/12). 
We model the dynamical noise by supposing that both 
6i{t) and Jij{t) change randomly after a time interval r. 
Within the time interval they are constant. 

For J ~ 6 ~ the spectrum of the Hamiltonian is 
composed of n-l- 1 degenerate levels, with interlevel spac- 
ing 2Ao, corresponding to the energy required to flip a 
single qubit. We study the case 6,J <^ Aq, in which 
the degeneracies are resolved and the spectrum is com- 
posed by n -|- 1 bands. In this limit the coupling between 
different bands is very weak. We assume free boundary 
conditions and express all the energy in units of Aq (we 
choose ?i = 1). 

In the following we analyze the behavior of the fi- 
delity [i| and error 

F{t) ^ mU{t)m^ , E = -\nF. (2) 

starting from an initial state |^) which is an eigenstate 
of CTp^ {j = 1, . . . , n), U{t) being the unitary evolution 
generated by (0. We concentrate on the central band 
of zero total magnetization which is characterized by the 
highest density of states and for which one expects that 
the effect of noise is most pronounced. 
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FIG. 1; 



Fidelity as a function of time for n 
-2 r — A in— 1\ 



14 qubits in 

the FGR regime {J = 2 ■ 10"^, 5 = 4 • 10"^) and from top 
to bottom T = 1,3,5,10,20,25 (static imperfections). Inset: 
Fidelity as a function of time in the ergodic {J = 5 = 2 ■ 10~^, 
dashed line), and in the FGR regime (full line). 



Results - The decay of the fidelity due to static imper- 
fections is displayed in the inset of Fig. The system 
is characterized by two distinct dynamical regimes 
depending on the critical value Jc ~ S/n: the Fermi 
Golden Rule (FGR) (J < Jc) and the ergodic regime 
(J > Jc) i^l^. The FGR is characterized by a Lorentzian 
local density of states with width r^*^^. The ergodic 
regime is reached when all the levels inside the band 
participate to the dynamics; the local density of states 
coincides with the density of states and has a Gaussian 
shape with variance r°'^s. In both regimes the decay of 
the fidelity is the Fourier transform of the local density 
of states |9j and follows an exponential and a Gaussian 
decay with characteristic decay times Ffgr oc J^/6 and 
Ferg oc respectively (see Inset Fig.^l 

In the case of dynamical imperfections, different 
regimes emerge as a function of the frequency l/r. Be- 
low a critical timescale Tc the different behavior due to 
the ergodic and FGR regimes cannot be resolved any- 
more. This can be clearly seen in Figs. I1I2I A smoother 
crossover appears at a higher frequency l/rp (Fig. ^ 
when the noise frequency become comparable with the 
single qubit natural frequency (~ Aq). The error Et{T) 
at (fixed) time t tends to vanish as r decreases. 

The explicit calculation of the error to order yields 

Et{T)=AJ^a^{Ng{T)+giAt)), (3) 
where t = Nt -f At, with N integer, < Ai < r, and 



gir) — 2 ds du sine {6u) [n^i -|- n|| cos(4Aou)] , 
Ja Jo 

(4) 

n^i {n^i) being the number of nearest-neighbor paral- 
lel (antiparallel) pairs in the initial state and sinc(a;) = 
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FIG. 2: Error _E as a function of r for t = 25, n = 10, J = 
5x 10~^, in the ergodic regime 5 = 5x10""^, n\i = 8 (squares), 
nil = 13 (triangles) and in the FGR regime 5 = 3 x 10~^ 
(circles). The fits are given by Eqs. and ll^J-® with cP' = 
1/12, ric = 13, Ao = 1 and (n||,n||) equal to (8,5) (dashed), 
(13, 0) (dot-dashed). The transition at is shown only in the 
former case. All the errors scale as (data not shown). 



sinx/x. The integration can be explicitly performed [UJ^ 
although the resulting analytic expression is not very 
transparent. Note that, due to the convexity of .g(r), the 
error Et{T) < AJ^a^tg{T)/T, being equal when t/r = N, 
thus providing a simple interpolation of (PJ. Moreover, 
the function g{T) can be approximated in several impor- 



tant limits. For tS <ti 1, sine (Su) ~ 1, whence 



9ir) 



i||Sinc^(2Aor)] 



(5) 



which yields 5(t) ~ tIcT^ 



7r/4Ao and gir) 



n^ir'^ (ergodic regime) for t > Tp (see Fig.EJ, where the 
total number of links Uc = n|| +?^||, unlike n|| and n^i, 
does not depend on the initial state j^*). On the other 
hand, when tS 3> 1, smc\Su) ~ jS{u) (FGR regime) 
and I0J reads 



nu^[Sr-H2ST) 



7 - 



(6) 



where 7 ~ 0.577 is Euler's constant. Substituting these 
approximate expressions in Eq. Q, the error at a fixed 
time t for different r values scales like 



Etir) ~ AJ^aH 




(all regimes) 
(all regimes) 
(ergodic) 
(FGR) 

(7) 

In Fig. 121 we show the scaling of Et{T) with r for different 
values of J. For the ergodic regime we choose J = S, 
while the FGR is characterized by J ^ (5. As r < Tc the 
two distinct ergodic and FGR behaviors of the static case 
(compared in Fig. |21only for the sets with n|| = 8) are 
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FIG. 3: Error at time t = 50, for n = 10, J = 5 ■ 10"^ and 
different 5 values. The squares represent the ergodic regime 
5 = J. The FGR regime is plotted for S = 1,2,3,5 • 10"^ 
(empty, pointed, dashed, full circles respectively). Inset: 
as a function of 6 for n — 10, 12, 14 (circles, squares and dia- 
monds respectively) . The dashed line is proportional to 6~^^^ , 
in agreement with Eq. @. 



not resolved. Equations © and ©, plotted in Fig.|21 are 
in excellent agreement with the numerical results. The 
additional kink at r ~ Tp = 7r/4Ao sets in when single 
spin dynamics starts to play a role. We also checked 
that Tp is independent on J and 6, in agreement with 
Eq. (jSJ. The transition at r = Tc is striking and occurs 
when the error starts deviating from the linear behavior 
given by Eq. {Tj). In fact, the crossover between the two 
regimes could be defined by equating the third and the 
fourth line of Q, that is for t — tt/6, which for S — 0.3 
would give r ~ 10.5. However, since the saturation value 
Et^r) = AJ'^a'^n-^iTTt/S given by Eq. ^ is reached only 
for 6t ^ 1 and since the transition is sharp, a much more 
accurate way to define Tc is by looking at the point for 
which the deviation from the linear behavior [third line 
in Eq. ^] becomes apparent. To this purpose we keep 
the next-leading correction to Eq. (O and approximate 
sinc^Su) ~ 1 - (5u)^/3 (for t <1/S) in the integral Q). 
For T > Tp we obtain 



{5rf - 



{5t 



i41 



18 



(8) 



If the plot resolution in Fig.[21is some fraction e of the to- 
tal vertical range 4J^cr^n||t^, the error curve starts devi- 
ating from the linear behavior when {t / t){5t)'^ / {1^5"^) ~ 



I.e. 



(18rt)! 



(9) 



which for t — 25, 5 ~ 0.3 and e — 1/40 yields Tc — 5, in 
full agreement with Fig. El In Fig. |21we show the error 
£'t(r) with fixed J and different 5 values. The scaling of 
the critical threshold Tc is clearly visible. We also checked 



that Tc does not depend on J (data not shown) . The inset 
of Fig. shows the dependence of Tc as a function of 5, 
confirming the prediction 

Theorem - After having presented the overall picture of 
dynamical imperfections on the fidelity of computation, 
we complete our analysis and set up a general framework 
to consider the effect of a time-dependent noise on the 
evolution of a quantum system 



(10) 



where is time independent and H{t) varies with a 
given characteristic time t, according to the stochas- 
tic process with independent increments ^(t) = 
X^fcLi X[fer-T.fer)(^) ^k, where xa is the characteristic 
function of the set A and {S^k\k are independent and 
identically distributed random variables, with expecta- 
tions E[£,k\ = 0, Var[^fc] = E[el] ^ < oo. The time 
evolution operator over the total time t = tN is given by 



t/r 



UN{t) = n exp[-i(i?o + ikV)T], 



(11) 



fe=i 



where a time-ordered product is understood, with earlier 
times (lower k) at the right. Let us assume, for sim- 
plicity, that and V are bounded operators, so that 
U{t) is a norm-continuous one-parameter group of uni- 
taries and all our subsequent estimates are valid in norm. 
We are interested in the existence and form of the lim- 
iting time evolution operator C/Ar(t) for iV ^ cx). When 
expanding the product, one finds that the term inde- 
pendent of i is 1, while the term proportional to t reads 
—iHot—iVt ik/N . Now, according to the weak law 

of large numbers [ill, f-limAr^oo SfeLi ^k/N = E[S,k] = 
0, for we assumed =0"^ < oo, and the limit is taken 

in probability. Therefore, for iV — > oo 



1 ^ 



(12) 



fc=i 



Analogously, by using the weak law of large numbers, one 
can prove that all higher powers of Vt vanish in the limit, 
thus obtaining 

U{t) = P- lim UN{t) = cxp{-iHot), (13) 

N^oo 



in the following sense 



N 



lim P i\\UN{t) - e^p{~iHot)\\ > e) = 0, 



(14) 



uniformly in each compact time interval. If the term 
£,{t)V is viewed as exemplifying the effect of (dynamical) 
error-inducing disturbances, the above result physically 
implies that the effects of the errors are wiped out if 
their characteristic frequency is sufficiently fast. This 
defines the purely dynamical regime. 



4 



Another viewpoint can also be adopted, that is some- 
what complementary to the above one. Given a charac- 
teristic frequency of the noise, it is possible to establish 
an effective value of the strength of the imperfections 
so that the above result holds (approximately). In this 
sense, a natural question is what happens for large but 
finite N. This (physical) question can be answered by 
remembering that under the same hypotheses, according 
to the central limit theorem, the limiting random vari- 
able 77 — limAr_^oo X]fc=i Cfc/v^ exists and is Gaussian 
with mean E[r]] = and variance E[ri'^] = cr^, namely it 
is distributed like / (77) = (27r(T^)^ 2 exp(— ry^/2cr^). Thus, 
by following the same steps that led to we find that 
for > 1 



(15) 



UN{t) - exp [-^ (iJo + vV/y/N^ t 



Equation (|15|l implies then that for fixed t, the sys- 
tem "feels" an effective interaction strength e^g = 
cr||T^||/V^cx cr||V^||V^. 

For intermediate values of N, Eq. (|15|l is no longer 
valid, because it hinges upon the commutativity of Hq 
and V. However, by assuming that V ^ Hq (e.g. in 
norm), a straightforward expansion shows that the per- 
turbation V is replaced by 



V{t) 



1 



iHat-ir-iHnt 



dt e^^^'We 



(16) 



so that, for T||iJo|| ^ 27r, the effective perturbation be- 
comes 



V{t) ^ Vz=Y.PkVPk. 



(17) 



where Pk are the eigenprojections of iJo {Hq = Yl, ^kPk)- 
This phenomenon is reminiscent of the quantum Zeno 
subspaces ^12] . 

The generalization of the above results to a Hamilto- 
nian with a family of independent stochastic processes 
with zero mean and finite variances is straightforward. 
This is the case of the Hamiltonian |[[J, which reads 

H{t)^Ho + S^oit)-Vo + 2jm-y, (18) 



where Hq = E, AoaP\ (Vo), = al", {V) 



J]) 



and = (^o)i and = = l,---,n) are 

independent random variables uniformly distributed in 
the interval [-1/2,1/2]. 

We can then reinterpret our previous results in the 
light of the above theorem, by applying the well-known 
static results to the (static) evolution with renormal- 
ized couplings (|15|l (with r/V t] ■ V). Thus, indepen- 
dently of the interaction strength and the correspondent 
dynamical regime, there is a quadratic decay law for suf- 
ficiently large N — t/r (or small r). 



Et{r) 



(r < Tp), 



(19) 



where t^=^ = AJ'^{^{r] ■ Vf\^) = AJ^Uca^ and ~ 
A(^^ [the i?o-tiniescale, see (|16|l ]. On the other hand, for 
smaller A^, i.e. r > Tj,, the effective interaction (I17f) is 



given by {Vz)ij = cr^'a 



^W^O')^ whence 



Etir) - ^rorgt^ = Tcrgir (t > Tp), 



(20) 



where Torg = 4j2(*|(r7 • Vz)^|*) = AJ'^m^ia'^ . There- 
fore, we recover the linear growth of the error (with the 
correct coefhcients), that describes both regimes up to 
in Eq. 0. 

Conclusions - We studied the effects of dynamical im- 
perfections on a general model of a quantum computer 
and identified several dynamical regimes, depending on 
the frequency of the external noise as compared with the 
coupling constants of the quantum computer. Above 
a threshold frequency, imperfections can be treated as 
static imperfections, although with renormalized param- 
eters. Below this threshold the different dynamical 
regimes induced by the presence of imperfections are not 
resolved. These results give a better comprehension of 
the general problem of noise in quantum computers and 
might suggest new strategies to develop general error cor- 
recting techniques. 
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